Biased tomography schemes: an objective approach 
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(Dated:) 

We report on an intrinsic relationship between the maximum-likelihood quantum-state estimation 
and the representation of the signal. A quantum analogy of the transfer function determines the 
space where the reconstruction should be done without the need for any ad hoc truncations of the 
Hilbert space. An illustration of this method is provided by a simple yet practically important 
tomography of an optical signal registered by realistic binary detectors. 
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The development of effective and robust methods of 
quantum state reconstruction is a task of crucial impor- 
tance for quantum optics and information. Such meth- 
ods are needed for quantum diagnostics: for the veri- 
fication of quantum state preparation, for the analysis 
of quantum dynamics and decoherence, and for informa- 
tion retrieval. Since the original proposal for quantum 
tomography and its experimental verification 0, this 
discipline has recorded significant progress and is con- 
sidered as a routine experimental technique nowadays. 
Reconstruction has been successfully applied to probing 
the structure of entangled states of light and ions, op- 
erations (quantum gates) with entangled states of light 
and ions or internal angular momentum structure of cor- 
related beams, just to mention a few examples 0. 

All these applications exhibit common features. Any 
successful quantum tomography scheme relies on three 
key ingredients: on the availability of a particular tomo- 
graphically complete measurement, on a suitable repre- 
sentation of quantum states, and on an adequate mathe- 
matical algorithm for inverting the measured data. In ad- 
dition, the entire reconstruction scheme must be robust 
with respect to noise. In real experiments the presence 
of noise is unavoidable due to losses and due to the fact 
that detectors are not ideal. The presence of losses poses 
a limit on the accuracy of a reconstruction. However, 
the very presence of losses can be turned into advantage 
and used for the reconstruction purposes. As has been 
predicted in Ref. p|, imperfect detectors, which are able 
to distinguish only between the presence and absence of 
signal (binary detectors) provide sufficient data for the 
reconstruction of the quantum state of a light mode pro- 
vided their quantum efficiencies are less than 100%. The 
presence of losses is thus a necessary condition for a suc- 
cessful reconstruction: An ideal binary detector would 
measure only the probability of finding the signal in the 
vacuum state. 

The required robustness of a tomography scheme with 
respect to noise is often difficult to meet especially if it is 
biased, that is, if some aspects of the quantum systems 
in question are observed more efficiently than the others. 
Since our ability to design and control measurements is 
severely limited, this situation will typically arise when 



one wants to characterize a system with a large num- 
ber or infinitely many degrees of freedom, for instance in 
the quantum tomography of light mode mentioned above. 
The standard approach is to truncate the Hilbert space 
by a certain cut-off, reducing drastically the number of 
parameters involved Needless to say, such ad hoc 
truncation lacks physical foundation. It may have bad 
impact on the accuracy of reconstruction or conversely it 
may lead to more regular results. The latter case may 
easily happen when an experimentalist seeks for the re- 
sult in the neighborhood of the true state. Such a tacitly 
accepted assumption may appear as crucial as it allows 
elimination of the infinite number of unwanted free pa- 
rameters. This drawback erodes the notion of tomogra- 
phy as an objective scheme. 

In this Letter we will propose a reconstruction proce- 
dure that is optimized with respect to the experimental 
set-up, representation and inversion, designed for dealing 
with biased tomography schemes. The recommended ap- 
proach to the generic problem of quantum state tomogra- 
phy will be demonstrated on the scheme of a light mode 
adopting elements of linear optics (beam splitter) with 
realistic binary detectors detecting the presence or ab- 
sence of the signal only. In addition, we will, for the first 
time present a statistically correct description of such a 
tomographic scheme. 

Let us develop a generic formalism for the maximum- 
likelihood (ML) inversion of the measured data. Let us 
assume detections of a signal enumerated by the generic 
index j. Their probabilities are predicted by Quan- 
tum theory by means of positive-operator-valued mea- 
sure (POVM) elements Aj, 



Pi = TrIAjP], 



< Aj < 1, 



(1) 



p being the quantum state. The observations Aj are as- 
sumed to be tomographically complete in the Hilbert sub- 
space we are interested in. No other specific assumptions 
about the operators Aj, their commutation relations or 
group properties will be made. In general, probabilities 
Pj are not normalized to one as the operator sum 



5> J= G>0 



(2) 
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may differ from the identity operator. Theoretical prob- 
abilities pj can be sampled experimentally by means of 
registered data Nj . The aim is to find the quantum state 
p from data Nj . 

The ML scenario hinges upon a likelihood functional 
associated with the statistics of the experiment. In the 
following, we will adopt the generic form of likelihood for 
un-normalized probabilities [6( 



log £ = 5^- log 



Pj 



(3) 



which should be maximized with respect to p. Here 
the index j runs over all registered data. The extremal 
equation for the maximum-likely state can be derived 
in three steps: (i) The positivity of p is made explicit 
by decomposing it as p = o^o. (ii) Likelihood © 
is varied with respect to independent matrix a using 
5(logpj)/5o~ — Aja^ /pj\ (iii) Obtained variation is set 
equal to zero and multiplied from right side by a with 
the result 



R P = Gp, R = Y1 



E/ Pj' Nj 



(4) 



where the operator G is defined by Eq. J5J and operator 
R depends on the particular choice of C. Notice that 
this equation may be cast in the form of Expectation- 
Maximization (EM) algorithm 



RgPg — Pg, 



(5) 



where R G = G~ 1 I 2 RG~ 1 / 2 and Pg = G 1 / 2 P G 1 / 2 . This 
extremal equation may be solved by iterations in a fixed 
orthogonal basis. Keeping the positive semi-defmitcness 
of pa [by combining Eq. ||3J with its Hermitian conjugate] 
the (n + l)th iteration reads 



Pg 



(n+l) 



= R^p^R^K R { g ] = G-^Rip^G- 1 ' 2 , 



? («) 



Starting with some initial guess Pq the iterations are 
repeated until the fixed point is reached. In terms of pa, 
the desired solution is then given by 



p = G- 1 ' 2 p G G 



-1/2 



(6) 



Going back to likelihood in Eq. J3J) we now see, that 
the operator G coming from the mutual normalization 
of probabilities, EjPj = Tr[/oG], provides a complete 
(normalized) POVM, which is equivalent to the origi- 
nal biased observations Af ^2,-G~ x l 2 AjG~ x l 2 = Iq- 
This establishes the preferred basis for a reconstruction. 
Due to the division by the operator G in Eq. © and in 
the sentence above the reconstruction can be done only 
in the subspace spanned by the non-zero eigenvalues of 
G. The spectrum of G plays therefore the role of to- 
mographic transfer function analogously to the transfer 
function in optical imaging. It quantifies the resolution 



of the reconstruction in the Hilbert space. Large eigen- 
value of G indicates that many observations overlapped 
in the corresponding Hilbert subspace and this part of 
the Hilbert space is more visible. The Hilbert subspace 
where the reconstruction was done is clearly not a subject 
of a free choice in the proper statistical analysis. This is 
the main result of this Letter. This also gives a clue how 
to approximate the solution in the infinite dimensional 
case simply by taking the subspace corresponding to the 
dominant eigenvalues. The result of reconstruction can 
be easily checked in the preferred basis afterwards. If the 
reconstructed state exhibits dominant contributions for 
the components with relatively small eigenvalues of G, 
the result cannot be trusted. 

The essence of the correct reconstruction inhere in the 
following recommended scenario: After collecting all data 
the optimal basis for reconstruction is identified as eigen- 
vectors of G operator. The truncation is achieved by tak- 
ing into account only those with dominant eigenvalues, 
where the ML extremal equation should be solved keep- 
ing the semi-positive definiteness of the density matrix. 
This establishes the quantum tomography as an objective 
tool for the analysis of infinite dimensional quantum sys- 
tems. Indeed, previously reported results of tomographic 
schemes have always considered the space for reconstruc- 
tion ad hoc: If one knows what the result should be it is 
not really difficult to get it. 

Let us illustrate this procedure on the following simple 
realistic detection set-up: the signal state (described by 
the density matrix p) of the input mode a is mixed on 
a beam-splitter with the probe coherent state \(3) of the 
mode b and the mixed field is detected on a single on/off 
detector. Then the probability p of having no counts on 
the detector is measured. 

Such non-ideal measurements have already been used 
for tomography purposes. The inference of a photon 
number distribution was proposed in [8( and experimen- 
tally realized in 0. A more advanced setup based on a 
multichannel fiber loop detector was developed and ex- 
perimentally verified earlier in |Toj . As proposed in [T]| 
and ^3 1 the reconstruction of a full density matrix can 
be done by measuring a coherently shifted signal. This 
reconstruction technique has also been implemented ex- 
perimentally as a direct counting of Wigner function |13| . 
However, the algorithms used for the quantum state re- 
construction were not robust as indicated by the fact that 
they could give non-physical results. This is due to the 
a priori constraints put on a quantum object, namely 
the semi-positive definiteness of a density matrix p > 0, 
which is not guaranteed in the above mentioned schemes. 
While it seems to be intractable to implement the con- 
dition of positive semi-definiteness in Wigner represen- 
tation, it can be done in the general formalism adopting 
the maximum-likelihood estimation. 

The probability of registering no counts on the detector 
is given by Mandel's formula 



p = (: cxp {-i/ c c f c} :), 



(7) 
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FIG. 1: Eigenvalues of the matrix G 10) truncated at 
N tr = 15. The simulated measurement was done at N p 
7-points equidistantly distributed in regions: (a) and (b) 
fle( 7 ) G [-2,2], 7m( 7 ) € [-2,2]; (c) Re( 7 ) G [-1,1], 
Im(rt) = 0; (d) i?e( 7 ) G [1, 1.01], Ira,( 7 ) = 0. In all panels, 10 
equidistant values of the detector efficiency were chosen from 
the interval r\ G [0.1,0.9]. 



where v c is the efficiency of the detector; and c are cre- 
ation and annihilation operators of the output mode, and 
:: denotes the normal ordering. For simplicity, we assume 
here that in the absence of the signal the detector does 
not produce any clicks; dark count are ignored. Let us 
assume, that the beam-splitter transforms input modes 
a and b in the following way: c = acos(a) + 6sin(o:). 
Averaging over the probe mode 6, from Eqs. Q one 
obtains 



p=Y J ^-v) n (n\D\ 1 )pD{ 1 )\n) 



(8) 



n=0 



where v — ^ c cos 2 (a), 7 = — /3tan(a), -D( 7 ) = 
exp { 7 a^ — 7*0} is the coherent shift operator, and |n) 
denotes a Fock state of the signal mode a. Using the 
operator notation R„. 7 = D(7)|n)(n|-Dl( 7 ), the no- 
count probability is generated by the POVM elements 
A„ :7 = — u) n 'R rh y and, defining a collective in- 

dex j = {v, 7}, the counted probability coincides with 
Eq. JTJ. 

Figure shows how a suitable choice of 7-points for a 
fixed truncation number Ntr = 15 can be achieved. Ob- 
viously, the amount of data used in Fig.^a) as compared 
to Fig. ^b) is excessive for the reconstruction. On the 
other hand, when the number of points is too small, or 
they are chosen in an inappropriate way, eigenvalues of G 
differ strongly making reconstruction unfeasible. For ex- 
ample, in Figure (d) the last eigenvalue is only ~ 10~ 5 . 
However, one needs to mention that the analysis of G 
provides a necessary but not sufficient condition of the 
reconstruction feasibility. In particular, a single 7 point 
measurement is not sufficient (just like a measurement 



FIG. 2: A reconstruction of the state l)10[l according to 
procedure ©. The following measurements were used: 
Re(j) = -0.2,-0.1,0,0.1,0.2; Im{i) = 0.1,-0.5,0,0.5,0.1; 
20 equidistantly distributed detector efficiencies in the inter- 
val [0.1,0.9] were used; the Hilbert space was truncated at 
Ntr = 5. Panel (a) shows the eigenvalues of the matrix G. 
Panels (b) and (d) show the real and imaginary parts of the 
reconstructed matrix (in Fock basis). They were obtained 
using 10 6 iteration of the EM algorithm. Panel (c) shows 
the variances of the real part (n < m) and imaginary part 
(n > m) of the reconstructed elements given by Eq. ©. 



in 7 = is able to give only the diagonal elements). 
One needs to measure in at least two different non-zero 
7 points. The confidence interval on the reconstructed 
density matrix elements can be provided with help of 
variance a(p mn ) = (F(p mn )N mes ) , where N mes is 
the total number of measurements, and the Fisher in- 
formation F can be defined for real part of the density 
matrix elements as 171: 



F(Re[p mn ]) = Y, 



d 



Pj 



dRe{p mn ) J2f Pj' 



(9) 



and similarly with Re changed to Im for imaginary part 
of p. 

To illustrate our discussion, let us consider a recon- 
struction of the following state (Figure EJ): 

|<£) = (|0)+exp{0.5i}|2))A/2. (10) 

The simulation was done using a total of 10 7 measure- 
ments collected in five different points on the phase plane 
7. In Fig. Ufa) one can see the eigenvalues of the matrix 
G (J2J. Obviously, the chosen set of points is suitable for 
the reconstruction. Notice the correlation between de- 
creasing eigenvalues and increasing errors in Figs. Eja) 
and (c). 

This objective approach may compared with alterna- 
tive schemes based on the reconstruction of Wigner func- 
tion. Measurement in any given 7 point is able to give a 
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FIG. 3: A reconstruction of the signal coherent state a = 
exp{j7r/4}; (a) the reconstructed Wigner function; (b) the di- 
agonal elements of the reconstructed density matrix; (c) the 
difference between the exact and the reconstructed Wigner 
functions; (d) the variance (7(7,7*). The Wigner function 
was reconstructed point- wise at 2500 points of the phase plane 
from N r = 10 4 measurements per point using Nu = 10 3 iter- 
ations of the EM algorithm. The Hilbert space was truncated 
at N tr = 12; 30 different values of detector efficiencies were 
used equidistantly distributed in the interval [0.1, 0.9]. 



value of the Wigner function in that point. Indeed |15| . 

W( 1 ) = -J2(-^) n Rn n , (11) 

where i?„( 7 ) = Tr[pR„, 7 ] = (n\D^(j)pD(j)\n). For a 



fixed value of the amplitude 7 one should seek the set of 
non-negative matrix elements i?„ i7 and plug in these val- 
ues into the definition of the Wigner function pi(l . These 
matrix elements can be found by inverting the counted 
statistics (JSJ measured with a set of different efficiencies 
solving a linear positive inverse problem. This can be 
accomplished by means of the EM algorithm similarly 
to the approach used in |lo|. An example of such a re- 
construction is shown in Fig. [3] Though the reconstruc- 
tion seems to be faithful, one should keep in mind that 
even very small deviations from the true Wigner func- 
tion might make it non-physical. Such Wigner function 
would not correspond to any physical, positive definite 
density matrix. This is due to the fact that the operators 
R n ,7 do not commute for different 7s, so noisy measure- 
ments may give inconsistent results. Going back from 
Wigner function to the density matrix using Glauber's 
formula [3, p = 2 J d 2 j{~l) n W(j* ,j)D(2j), one can 
see in Fig. E^b) , that some diagonal elements of the re- 
constructed matrix are negative. 

A generic biased tomography scheme addressing some 
aspects of the quantum systems more efficiently than 
other aspects was introduced. Its performance is charac- 
terized by quantum analogy of transfer function, which 
may be further optimized to achieve the desired resolu- 
tion. This establishes tomography as an objective tool for 
quantum diagnostics. The recommended approach was 
demonstrated on a simple, robust and effective quantum 
tomography scheme using detectors that are only capa- 
ble to distinguish between the presence and absence of 
photons. 
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